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In the diffusion-collision model, the unfolding rates are given by the likelihood of secondary struc- 
tural cluster dissociation. In this work, we introduce an unfolding rate calculation for proteins 
whose secondary structural elements are a-helices, modeled from thermal escape over a barrier 
which arises from the free energy in buried hydrophobic residues. Our results are in good agreement 
with currently accepted values for the attempt rate. 

PACS number(s): 87.14.Ee, 87.15. Aa, 82.20.Db 



I. INTRODUCTION 

In the diffusion-collision model of protein folding ^ the 
protein is modeled using a collection of spheres connected 
by flexible strings. The spheres represent the secondary 
structural elements such as a-helices or /3-sheets (or clus- 
ters of these secondary structures) called microdomains, 
that constitute the protein. 

The folding process from a completely unfolded protein 
to the the final native state is accomplished via diffusion 
through the solvent, collision, and finally coalescence of 
the microdomains. The state of the protein is defined by 
the number of pairings between the microdomains that 
are present at a given time t. The rate equations for 
transitions between these states can be written as 



dPjt) 
dt 



= KP{t) 



(1) 



where P(t) is the vector of states and K is a. matrix con- 
taining the transition rates between the different states. 
A protein having, say, q microdomains would involve 
p ^ q{q - l)/2 pairings, 2^ states P^{t) and a 2'' x 2^ 
rate matrix K. 

In general, the calculation of the elements of the rate 
matrix K is somewhat involved. The forward rates are 
the rates of structural coalescence. In the diffusion- 
collision model the forward rates are calculated assum- 
ing the microdomains diffuse through a solvent environ- 
ment, the space of which is limited by the length of the 
intervening strings and the van der Waals radii of the 
microdomains. These microdomains are assumed to be 
nascently formed, and their degree of formation is given 
by a helix-coil transition theory calculation |^ (as in 
AGADIR [HQl) in the case of a-helices, or via a combi- 
nation of theory Q and experiment Q in the case of /3- 
sheets. As the microdomains undergo diffusion, they oc- 
casionally collide. When this happens the microdomains 
coalesce with a probability 7, being held together by hy- 
drophobic interactions in the case of a-helices, or a com- 
bination of hydrophobic and hydrogen bond interactions 
in the case of /3-sheets. The coalescence probability 7 
is given by the likelihood that the microdomain is in a- 
helical or /3-sheet form, the percentage of hydrophobic 



area, and the likelihood of proper geometrical orienta- 
tion upon collision. 

The forward folding times between any two given states 
in the mean first passage time approximation M are given 
by 



(2) 



where V is the diffusion volume available to the mi- 
crodomain pair, A is the target area for collisions, D 
is the relative diffusion coefficient, 7 is the probability of 
coalescence upon collision and I and L are geometrical 
parameters calculated for diffusion in a spherical space. 
The inverse of the first passage time-scales t/ are the 
forward folding rates kf that are used in the rate matrix 
K. 

The pairs can also dissociate. In typical diffusion- 
collision model calculations, the form of the unfolding 
times T(, used for two microdomains A and B comes from 
the Van't Hoff-Arrhenius law given by 



(3) 



where AGab is the free energy difference between paired 
and unpaired states, ks is Boltzmann's constant, T is 
the temperature and v is an attempt rate. In the case 
of a-helices the dominant contribution to the free energy 
comes from the buried hydrophobic area and therefore 



AGab — IAab 



(4) 



where / is the free energy change per unit buried hy- 
drophobic area in the pairing |^ and Aab is the buried 
area The unfolding rates fcf, are given by the inverse 
of the unfolding times Tjy. 

The diffusion-collision model has been successful in de- 
scribing the overall folding kinetics of several proteins 
|]lO|-p^ . In each of these studies a single value of the pa- 
rameter V was used for every unfolding transition. This 
value was adjusted to obtain the desired result, namely, 
to ensure that the protein would fold to its native state. 
This procedure is justified because the equilibrium (or na- 
tive) occupation probabilities are known; in fact, for suffi- 
ciently simple systems the folding and unfolding rates can 
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be determined from these probabilities ||T^. The typical 
values used lie between 1-1000 ns~^, which yields unfold- 
ing rates consistent with observed rates of bi-molecular 
dissociation [jl^ . 

In cases where the final occupation probabilities are 
unknown, for instance in studies of protein mis-folding 
and non-native kinetic intermediates [T5| ] such methods 
are clearly not possible. Indeed, even a detailed de- 
scription of the intermediate folding kinetics of a protein 
whose final state is known requires a more accurate and 
foundational determination of u, as was pointed out by 
Burton et al. 

In this work we compute unfolding rates that, in the 
context of the diffusion-collision model, can be used for 
any given unfolding transition in the study of proteins 
whose secondary structural elements are a-helices. From 
the rates we find the values of the parameter v. This 
makes the diffusion-collision model more predictive and 
enables it to be used in situations where the occupation 
probabilities are unknown. 




FIG. 1. Potential for the two microdomains A and B. The 
potential is infinite on the left because of the hard-core repul- 
sion of the van der Walls contact between the microdomains. 
The barrier on the right can be crossed by microdomain pairs 
with energies larger than Et = /Aab, the free energy difi'er- 
ence between paired and unpaired states with a buried hy- 
drophobic area Aab- The width of the well L, is taken to be 
the diameter of a water molecule. 



II. CALCULATION OF THE UNFOLDING RATES 

We model the dissociation of microdomains as a ther- 
mal escape event over a barrier. Consider the pair of 
microdomains (which could be a-helices or clusters of a- 
helices) A and B connected by a string, diffusing in the 
potential well depicted in Figure |l|. The left boundary 
is infinite because of the hard-core repulsion of the van 
der Waals contact between the pair. Pairs with energies 
larger than Ei, = /Aab, the free energy difference be- 
tween paired and unpaired states, can escape from the 
right boundary of the well. The well width L is set to 
the diameter of a water molecule. A separation larger 
than L exposes the buried hydrophobic area of the pair 
to the solvent, the free energy savings is lost, and the 
pair separates, resulting in an escape from the potential 
well. 

The binding energies Ei, of microdomain pairs in pro- 
teins are typically much larger than the thermal energy 



Eb > keT. 



(5) 



This means that the time to escape from the well is much 
larger than any other time-scale involved in the problem, 
in particular larger than the thermalisation (or velocity 
auto-correlation) time and larger than the time it takes 
for the pair to diffuse in the well. Consequently, at any 
one time, the spatial distribution inside the well of an 
ensemble of pairs will be homogeneous 



p{x, t) (X 1/L 



(6) 



and the flux incident on the barriers will be thermal. We 
will use these two facts to calculate the rate at which the 
pairs dissociate. 



The flux at the boundary on the right (at x = L) 
depends on the density of pairs at that boundary and the 
probability that their energy is high enough to thermally 
escape over the boundary. The differential element of flux 
at the boundary L of pairs with relative velocity between 
V and V + dv is given by 



dJ°''\L,t) ^ vp{L,t)dN{v) 



(7) 



where p(L, t) is the number density of pairs at the bound- 
ary at a time i. 



dN{v) 



27rkBT 



1/2 



(8) 



is the fraction of pairs with relative velocities between v 
and V + dv, and fi is the reduced mass given by 



TTIA + TUB 



(9) 



where ruA and tob are the masses of the two mi- 
crodomains. 

In order to find the total fiux through the outer bound- 
ary at L we must integrate over all velocities larger than 
+ -\/i?b/2r7i since the potential barrier can be crossed by 
pairs with energies higher than _Et,, and pairs with rela- 
tive velocities higher than that can contribute to the flux 
leaving the well. This yields a flux out of the well 



r^\L,t) = p{L,t) 



(hT\ 

\27TpJ 



1/2 



-Et/keT 



(10) 



If the number of pairs inside the well at some time t is 
n{t) then, because of (|^), the number density must be 
p{x,t) = n{t)/L everywhere and 
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n{t) fhT\ 

L \2TTflJ 



1/2 



(11) 



This means that the dissociation rate constant for a pair 
of microdomains with reduced mass /i and buried hy- 
drophobic area Aab = Ei, / f at a temperature T is 



1/2 



(12) 



The terms preceding the exponential correspond to our 
prediction for the Van't HofT-Arrhenius attempt rate v 
in (^). As an example, the attempt rate found for a 
coalesced pair of 16-residue Regan-Degrado [M helices 
with a combined hydrophobic area loss of 600A^ is 64 x 
lO^s-i. 

It is interesting to note that a result similar to ([l^ ) 
would have been obtained by assuming the the attempt 
rate to be the inverse of the thermal well-crossing time, 
namely taking 



1 knT 



(13) 



in ( |3|). This is, in fact, not the origin of the pre-factors 
in (pl2[). They arise as a consequence of (^): The fac- 
tor of 1/L comes from the homogeneity of the spatial 
distribution (^ and the factor of ^JksT /2iiii from the 
integration of the thermal velocity distribution (|^). 

It is possible that dissociation events within a pro- 
tein also include a relative rolling and/or sliding mo- 
tion of the microdomains. In this case the calculation 
above can be performed with a few minor differences 
that take into account the extra degrees of freedom. The 
relative velocity distribution of the microdomains is still 
the one-dimensional Maxwell-Boltzmann distribution be- 
cause motion parallel to the surface through which the 
probability is flowing does not contribute to escape from 
the well. The probability in the bound region is homoge- 
neously distributed in a two- or three-dimensional volume 
in ( p^ ) , and flows out of that volume through a one- or 
two-dimensional area. This calculation yields the result 



1/2 



(14) 



where we set d = 2 if we include either the rolling or 
sliding degrees of freedom, and d = 3 if both of them are 
included. Due to the steric clashing of the side chains 
it seems rather unlikely that dissociation would include 
a sliding motion along the axes of the microdomains. It 
may be relevant, however, in the context of molten glob- 
ules. 

This approach succeeds in removing the free parame- 
ter z/ from the diffusion-collision model. Moreover, our 
results for the one-, two- and three-dimensional unfold- 
ing rates have a \JT j dependence that could be used 
to distinguish between this and other proposals for the 
mechanism of microdomain pair dissociation. 



III. CONCLUDING REMARKS 

We have presented a calculation for the dissociation 
rate of a microdomain pair using a simple potential bar- 
rier over which pairs having energies above the free en- 
ergy of the hydrophobic docking can escape. Since we 
have not accounted for the energy in hydrogen bonds 
this result is relevant for the dissociation of a-helix pairs 
or helix cluster pairs only and not for the dissociation of 
/3-sheet pairs. We have found the unfolding rates arising 
from thermal fluctuations out of this potential well to be 
in good agreement with currently accepted values of the 
attempt rate v. 

The motivation of this work was to eliminate the free 
parameter v from the diffusion-collision model. In previ- 
ous applications of the diffusion-collision model (see for 
example JTo|-[l^ ) the folding kinetics from a denatured or 
random coil state to the final native state were followed. 
In such a case, it is reasonable to set the parameter v 
such that the native state achieves most of the proba- 
bility, because we know that the final state is attained 
at the end of the folding process. However, the removal 
of this parameter is important when considering folding 
processes which do not involve the native state. For ex- 
ample, in studying intermediate processes or protein mis- 
folding |l^, where the occupation probabilities may be 
completely unknown, such reasonable estimates of v are 
not available. In these cases, elimination of z^ as a free 
parameter is crucial. 

The results presented here also predict a oc \JT j ^ 
dependence in all cases which can be distinguished exper- 
imentally from other possibilities such as v ccT jl^ . An- 
other difference is the dependence of the unfolding rates 
on the states, not only through the hydrophobic area, 
but also through the reduced mass /i of the microdomains 
undergoing dissociation. This is markedly different from 
typical diffusion-collision model calculations where the 
attempt rate v is assumed to be the same for all dissoci- 
ation events within the protein. 
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